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E-8220 


ON THE KERNEL FUNCTION FOR THE UNSTEADY SUPERSONIC CASCADE 


WITH SUBSONIC LEADING EDGE LOCUS 
by M. E. Goldstein 
INTRODUCTION 

In current aircraft engine technology there has been considerable 
interest in the problem of the unsteady supersonic cascade with subsonic 
axial velocity* Thus we consider a two-dimensional oscillating cascade 
with a subsonic leading edge locus in a supersonic flow which is uniform 
far upstream. We suppose that the blades have small thickness and 
camber and are undergoing small amplitude harmonic oscillations. Kurosaka^ 
has obtained a low frequency analytical solution to this problem and Verdon^ 
has obtained a finite difference solution. 

In this note we reduce the problem to the solution of a functional 
integral equation and derive and compare two representations of the kernel 
function which are useful for computations, 

DERIVATION OF INTEGRAL EQUATION 

Consider a two-dimensional oscillating cascade with a subsonic lead- 
ing edge locus in a supersonic flow which is uniform far upstream. The 
blades which are assumed to have small thickness and camber, are undergo- 
ing small amplitude harmonic oscillations. We suppose that all lengths 
are nondimensionalized by the half-blade chord, c/2, the time t is non- 
dimensionalized with respect to c/2 divided by the free stream velocity 

U, the pressure fluctuation p is nondimensionalized by the free -stream 

2 

density times U , the upwash velocity v is nondimensionalized by U 

and a denotes the free stream speed of sound. Then the pressure fluctua- 
tion is governed by the equation (see fig. 1) 
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M » U/a is the free stream Mach number and the upwash velocity is 
related to the pressure by 
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The upwash velocity on. the n L blade is assumed to differ from that 
of the 0 blade by only a phase factor so that 


V(x + ns^ # ns) = e ina V(x,0) 


for x < 1, n » 0,±l,+2, ... (4) 


where a is the interblade phase angle and the upwash velocity on the 
0 th blade is related to its displacement Wge la3t by 


V(x 


,0) = (-i« + w 0 
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for x < 1 


As is usual we suppose for convenience that the frequency has a small 
positive imaginary part which we shall set equal to zero at the end of the 
analysis . Then 

k = k r 4- ie with 0 < e « 1 

and the outgoing wave boundary condition at infinity is now replaced by a 
boundedness condition. 

Since equation (1) possesses the separation of variables solution 
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the boundary condition (4) suggests that (following Lane and Friedman ) 
we seek a solution in the form of the superposition 
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we have put (see fig. 1) 

t 

v 

n = 0 , ±1 , ±2 , . . . 


x - x ~ ns 
n 


y = y - ns 

3 n 3 

and in order to insure that the solution remain bounded at infinity we have 
chosen the branch cut for the square root, Y, and the integration contour, 
a r +16, in the manner shown in figure 2 (with 6 > e). This solution 
possesses the jump discontinuity 

^\°°+i<$ 

[?(x)] - [<r n (x)] = ( f n (a)e' laXn da (7) 

v...'_oo+i(5 

across the line y = ns passing through the n th blade since it is only 
possible to satisfy the requirement that the upwash velocity be continuous 
by allowing a discontinuity in the pressure. The resulting pressure dis- 
continuity (in front of and behind the blade will be eliminated in the 
subsequent analysis. 



Since the upwash velocity, v, vanishes at infinity, equation (3) can 


be integrated to obtain 
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Inserting equations (5) and (6) and carrying out the integration now shows 
that 
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If we put 
£ n (a) = e inr f 0 (a) 
where 

T = a - Mks + 

It is easy to show from equation (8) that 
V(x + ns + , y + ns) = e ina V(x,y) 

Hence the boundary condition (4) is automatically satisfied and 
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On the other hand since [^(x)] = 0 for ]x| > 1 we can invert the 
Fourier transform in equation (7) (with n = 0) to obtain 
1 
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Upon inserting this into equation (10) and interchanging the order 


of integration, we obtain 

V(x,y) = K(x - x’ ,y) [P(x' ) ]dx' 

where 

p = pe^ = t e iMKX 
and 
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By letting y -> 0 we obtain an integral equation for the pressure 
ium [P] across the zeroth blade in terms of the known upwash velocity on 
the blade surface. Namely, 


V(x,0) - 
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where 

K (x) = lim K(x,y) (14) 

yH) 

EXPRESSIONS FOR THE KERNEL FUNCTION 

The form (12) and (14) for the kernel function is not suitable for 
numerical evaluation because the integral will not converge if we just put 
6 * e = 0 in the integrand. In order to carry out this limit It is con- 
venient to express the kernel in a different form. 

Form 1 . - Using the results of appendix A shows that 
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Because of the Heaviside function, H, the sum only need to be carried out until 

+ i 

the largest n for which ns T + £s|n| < x. For a subsonic leading edge s > s0 

i 

hence for |xj <1 no terms with n greater than l/(s + 0s) can occur 
in the sum. 

Form 2. - Since .^(a - Mk) * 0 and J/Ly > 0 for <5 - Me 
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and we can use the geometric series 
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where 
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At first glance it might appear that the integrand in this expres- 
sion possesses branch points due to the appearance of the radical y. 

However, it can easily be verified by replacing y by -y that this 

2 

function depends only on y and the branch points are therefore 
"canceled, 11 and the integrand possesses only poles. 

We can therefore use Jordan* s lemma to evaluate the integral in 
terms of its residues. To this end notice that the poles of the integrand 
occur at a * k/M and at the points where 
A + = nT r for n = 0, ±1, ±2 

But it follows from equation (16) that the latter points are given by 


( 18 ) 



where we have put 
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and 

d + = \/s +2 - g 2 s 2 (20) 

Notice that d is real for the subsonic leading edge. The + sign 
corresponds to the roots which lie in the upper half plane and the - sign 
to those in the lower half plane. The locus of roots in the complex 
a-plane is shown in figure 3. 

When x < 0 we must close the contour in the upper half plane and 
when x >0 in the lower half plane. Hence 
}K + (x) x < 0 

Kq(x) =( (21) 

1k”(x) x > 0 

where 


K^Cx) = ± Res in J half plane 


Then upon using the results of appendix B to evaluate the residues we find 
that 

( 22 ) 
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These series are only conditionally convergent and will not converge at 
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all if we take the derivative. In order to obtain convergent series 
notice that 

r 2 
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Hence the n t ^ 1 term of these sums behaves like 
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The series composed of these terms will converge to a row of step 


functions. Hence its derivatives will converge to a row of delta func- 


tion. We can evaluate the latter series by using the theory of distri- 
butions to show that (Lighthill^) 
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are now convergent series. 

COMPARISON OF KERNEL FUNCTION REPRESENTATIONS 

When the kernel function given by equations (21), (25), and (26) or 
that given by equation (15) is substituted into the integral equation (13) 
we obtain a functional integral equation (and not an ordinary integral 
equation) due to the introduction of terms of the form 
[p(x n + |n| s$) ] 

caused by the integration of the delta functions . The second form of the 
kernel function has two advantages over the first. Namely its series 
converge like e' Lnx /n rather than like e 1_Tlx //n and it is much simpler 
and faster to evaluate the terms of its series than evaluating Bessel 


functions . 
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The series in (25) and (26) are only conditionally convergent. 
However, it is shown in appendix C that the same device which was used 
to make the original converge can also be used to replace the series 
(25) and (26) by absolutely convergent series. Although the results, 
which are given by equations (Cl) and (C2) are now more complicated, 
they are definitely more suitable for numerical computation. They show 
that the removal of the slowly convergent part of the series results in 
a row of step functions which represent the discontinuities of K~. The 
series which appear in (Cl) and (C2) represent continuous function. 
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APPENDIX A 
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APPENDIX B 
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it follows that the residue of the integrand at is 
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APPENDIX C 

—1 

Upon retaining terms of order n we find that the n term of 
the sums in equations (22) and (23) now behaves like 

f + - 2n,>T' 0) [ s + * ^ )]} (X ? Py) ] 33 n * " 


Hence we can improve the convergence of the sums in equations (25) 
and (26) by adding and subtracting 
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Fig. 1 - Dimensionless Cascade Configuration. 





